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Abstract. We used four known chromospheric activity indicators to measure long- 
term activity variations in a sample of 23 M-dwarf stars from the HARPS planet search 
program. We compared the indices using weighted Pearson correlation coefficients 
and found that in general (i) the correlation between S can and Na i is very strong and 
does not depend on the activity level of the stars, (ii) the correlation between our S can 
and Ha seems to depend on the activity level of the stars, and (iii) there is no strong 
correlation between Scan and He i for these type of stars. 



1. Stellar magnetic cycles of M-dwarf stars 

There is currently a focus on the search for planets orbiting around M-dwarf stars. Due 
to their low mass, it will be easier to find lower mass planets orbiting these stars with 
the radial-velocity (RV) technique. And therefore it is extremely important to access 
all sources of intrinsic noise that can degrade the quality of the detected RV signals. 
There are hints that the magnetic cycles of stars may ind uce RV capable of hiding the 
signals of extra-solar planets (e.g. lDumusque et"aT1l2010h but not much is known about 



the long-term activity variations of M-dwarf stars. 

The Mt. Wilson surv ey showed that many solar-like star s have magnetic cycles 
similar to that of the Sun dWilson|[l978l : iBaliunas et al.| [l~995) but only one M dwarf 



was actually included in the sample. This star, HD 95735, was shown to have long- 
term activity variations without a defined cycle. More recent studies however uncovered 
evidence for the existence of periodicity in the long-term activity of a few M stars. Our 
closest neighbour, Pro x Centauri (dMe 5.5), was found to have a magnetic cycle with 
a -440 days period dCincunegui et al.ll2007al ). iDfaz et all d2007bh also found a -760 



days periodici ty in the activit y of the spectroscopic binary G1375 (dMe 3.5). More 



recently, iBuccino et al.1 d2010h announced the detection of cycles with periods of ~4 
and ~1 years for G1299A (Ml/2) and G1752A (M2.5), respectively. In a recent paper, 
we studied the influence of the long-term activity cycles in the RV signals of a s ample 
of 7 early-G and early-K stars known to have activity cycles (San tos et alj|2010l) . We 
found no hints of RV induced variations by the activity cycles of these stars at the ~lm 
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s level achieved by HARPS. In the present work we extend this study to the lower end 
of the main sequence by first analyzing the long-term behavior of the chromospheric 
activity indices and posteriorly compare them with the RV and parameters of the cross- 
correlation function of the stars. 



2. Sample 

We u sed a sample of stars from the HARPS M dwarf planet search program (iBonfils et al 



2010l) . This corresponds to a volume limited selection of stars closer than 1 1 pc, with a 



declination 8 < +20 degrees, brighter than V = 14 mag and with a projected rotational 
velocity vsin/ < 6.5 kms" 1 . The data consists of high resolution spectra which span 
around 6 years, from 2003 to 2009. 

Our first selection was to remove all spectra with S/N < 2 at the spectral order 6 in 
the region of the Ca n K line. All data was then nightly averaged and stars with more 
than 15 measurements were selected. We made bins of 150 days to average out possible 
rotationally modulated activity variations. For the rest of the study only stars with more 
than four bins were considered. This process resulted in a sample of 23 stars ranging in 
spectral type from M0 V to M4 V. 



3. Measuring the activity indices 

To measure the activity variations of the stars we computed four known chromospheric 
activity proxies based on the Ca n H and K, Ha, Na i Dl and D2, and He i D3 lines . 
All these indices were c alculated similarly to the procedure used in lSantos et all (120001) 
and lBoisse etail (l2009h . 



Our S call index was computed by measuring the flux in the core of the H (/13968.47) 
and K (/13933.66) lines using bands of 0.6 A and dividing it by two reference bands of 
20 A centered at 3900 and 4000 A. This index was then weighted by the square of the 
errors in each line taken as where N is the number of counts inside each band. The 
Ha index was computed in a similar fashion, by measuring the flux in a 1.6 A band cen- 
tered at the line core (/16562.808) and dividing it by two reference windows of 10.75 and 
8.75 A centered at 6550.87 and 6580.31 A, respectively. We also used a i ndex based 
on the Na i /15895.92 and /15889.95 lines similar to the one proposed by iDf az et all 
(2007a|). The flux in the cores of the two lines was measured using 0.5 A bands and 



divided by two reference regions of width 10 and 20 A centered at 5805.0 and 6090.0 
A, respectively. Finally, the He i index was obtained by integrating the flux in a 0.4 
A band at the line core at 5875.62 A and dividing it by 5 A reference bands at 5869.0 



and 5881.0 A as in lBoisse et all (I2009ri . 

The errors in the four indices were estimated by differentiating the respective equa- 
tions and taking into account the flux in each band used. 



4. Correlations between the activity indices 

We used the weighted Pearson correlation to compare our four different indices. All 
coefficients were calculated using the binned data. To access the significance of the ob- 
tained coefficients we computed the false-alarm-probability (FAP) by bootstraping the 
nightly averaged measurements, then binning the data and recalculating the correlation 
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coefficient for each of the 10 000 permutations. We chose significant FAPs those values 
smaller than 0.05, corresponding to a 95% significance level. The results are presented 
in Tabled] with the significant FAPs marked in bold. 



Table 1. Weighted Pearson correlation between S can and the other activity in- 
dices. FAPs calculated using bootstrap permutations (see text). On bold are FAP 
values lower than 0.05. 



Star 


p(5 Cfl //,Ha) 


FAP 


p(5 Cfl //,Na I) 


FAP 


p(Sc a ii, He I) 


FAP 


GJ361 


0.127 


0.4326 


0.744 


0.1335 


-0.086 


0.4686 


Gil 


-0.048 


0.3996 


0.965 


0.0076 


0.611 


0.1881 


Gil 76 


0.932 


0.0036 


0.413 


0.2206 


0.804 


0.0211 


G1205 


0.946 


0.0063 


0.861 


0.0167 


0.267 


0.3089 


G1273 


-0.119 


0.4432 


0.905 


0.0236 


0.386 


0.2757 


G1382 


0.997 


0.0001 


0.987 


0.0012 


0.679 


0.1162 


G1433 


-0.407 


0.2505 


0.985 


0.0006 


-0.370 


0.3275 


G1436 


-0.138 


0.3881 


0.860 


0.0124 


0.436 


0.1706 


G1479 


0.992 


0.0040 


0.945 


0.0245 


0.880 


0.0633 


G1526 


-0.823 


0.1031 


0.992 


0.0104 


0.243 


0.3520 


G1581 


-0.656 


0.0176 


0.653 


0.0269 


0.593 


0.0624 


G1588 


0.672 


0.1179 


0.953 


0.0173 


0.820 


0.0990 


G1667C 


-0.512 


0.0829 


0.863 


0.0012 


-0.487 


0.0660 


G1674 


0.411 


0.3113 


0.893 


0.0493 


0.640 


0.1998 


G1680 


0.919 


0.0383 


0.643 


0.1461 


0.754 


0.1492 


G1699 


0.035 


0.5079 


0.675 


0.2239 


-0.334 


0.3521 


G1832 


0.519 


0.1411 


0.976 


0.0003 


0.251 


0.3430 


G1849 


0.695 


0.0951 


0.723 


0.0557 


0.209 


0.3733 


G1876 


0.843 


0.0984 


1.000 


< 1.10 5 


0.771 


0.1343 


G1887 


0.831 


0.1112 


0.896 


0.0596 


0.571 


0.2291 


G1908 


-0.512 


0.1272 


0.927 


0.0056 


0.331 


0.2973 


HIP 12961 


0.383 


0.2153 


0.092 


0.4436 


0.700 


0.0738 


HIP85647 


0.955 


0.0112 


0.987 


0.0010 


0.914 


0.0235 



4.1. Scan versus Ha 

From our sample of 23 stars, 7 have significant correlations between the S call and the 
Ha index, representing ~30 % of the total (Table [lj. These stars are G1176, G1205, 
G1382, G1479, G1581, G1680, and HIP85647. O ne of them, G1581 prese nts anti- 
correlation between these indices. As was found by ICincunegui et all ( 2007bh for their 



study of 109 stars ranging from F6 to M5, we found a great range of correlations with 
coefficients from p = -0.823 (G1526) to p = 0.997 (G1382). Figure Q] shows what 
seems to be a trend between the correlation coefficient and the average Scan for each 
star. Stars with smaller average activity level tend to have negative correlations while 
for values of S call > 0.035 the correlations are all positive. 

4.2. S call versus Na i 

Figure |2] shows the weighted Pearson correlation coefficient between Scan and Na i 
versus the mean activity level measured by the Ca ii lines. It is notorious the strong 
correlations that can be observed between these two indices. There are no negative 
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Figure 1. Weighted Pearson correlation coefficient between Scan and Ha versus 
mean activity level. Black squares with FAP < 0.01, black points with 0.01 < FAP < 
0.05, small open circles are correlation coefficients with FAP > 0.05. 



coefficients and all but two stars (G1176 and HIP12961) have coefficients with values 
higher than 0.5. Furthermore, ~70% of the stars show correlations with FAPs lower 
than 0.05 and there seems to be no trend between the correlations and the average S call 
activity level. The stars with significative correlation coefficients (FAP < 0.05) are Gil, 
G1205, G1273, G1382, G1433, G1436, G1479, G1526, G1581, G1588, G1667C, G1674, 
G1832, G1876, G1908, and HIP85647. 

4.3. S can versus He i 

The correlation between these two indices is very week. As we can observe in Figure 
[3] there is a tendency for positive correlations as only four stars present negative coeffi- 
cients. But only two stars (G1176 and HIP85647, representing ~9% of the sample) have 
significant correlations between the two indices (with FAP < 0.05). Furthermore, there 
is a large dispersion of values of the correlation coefficient ranging from p = -0.487 
(G1667C) top = 0.914 (HIP85647). 



5. Conclusions 

We measured four activity indices for a sample of 23 M stars from the HARPS planet 
search program during a timespan of around 6 years. We compared the activity indices 
using weighted Pearson correlation coefficients and found that: 



There is a strong correlation between our S call and the Na i indices. This con- 
fir ms that the Na i line s are good activity proxies for these cool stars as suggested 
bv lDiazetal.1 (l2007ah for very active stars. 
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Figure 2. Weighted Pearson correlation coefficient between S call and Na i versus 
mean activity level. Symbols as in Fig. Q] 
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Figure 3. Weighted Pearson correlation coefficient between S can and He i versus 
mean activity level. Symbols as in Fig. Q] 



As observed by ICincunegui et all ( 2007bh we found a great range of correlations 
between the S can an d Ha indices. Furthermore we found what appears to be a 
trend between the correlation and the average activity level of the stars as mea- 
sured by the 5 call index. 
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Although some authors suggest the use of the He i line as a chromospheric ac- 



tivity proxy (e.g. ISaar et al.ll 19971) we found that this index is not well correlated 
with Scan for M dwarfs. 

Since the signal-to-noise ratio in the Ca u H & K lines is very low for M dwarfs, 
we suggest the use of the Na i Dl and D2 lines, situated in a redder spectral region, as 
an alternative chromospheric indicator for this type of stars. These results may influ- 
ence the way chromospheric activity is accessed in M-dwarf stars and contribute to the 
knowledge about the activity cycles of such stars. 

A more detailed study about this subject will be described in a future publication. 
Those results will then be used in the context of planet detection to search for trends 
between the long-term magnetic activity, RV, and parameters of the cross-correlation 
function in order to access at which level activity cycles might be inducing RV varia- 
tions for these type of stars. 
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